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This paper reports on a theoretical analysis of the rich variety of spatio-temporal pat¬ 
terns observed recently in inclined layer convection at medium Prandtl number when 
varying the inclination angle 7 and the Rayleigh number R. The present numerical in¬ 
vestigation of the inclined layer convection system is based on the standard Oberbeck- 
Boussinesq equations. The patterns are shown to originate from a complicated competi¬ 
tion of buoyancy-driven and shear-flow driven pattern forming mechanisms. The former 
are expressed as longitudinal convection rolls with their axes oriented parallel to the 
incline, the latter as perpendicular transverse rolls. Along with conventional methods to 
study roll patterns and their stability, we employ direct numerical simulations in large 
spatial domains, comparable with the experimental ones. As a result, we determine the 
phase diagram of the characteristic complex 3D convection patterns above onset of con¬ 
vection in the 7 — i? plane, and find that it compares very well with the experiments. In 
particular we demonstrate that interactions of specific Fourier modes, characterized by 
a resonant interaction of their wavevectors in the layer plane, are key to understanding 
the pattern morphologies. 


1. Introduction 

Pattern forming instabilities in macroscopic dissipative systems, driven out of equilib¬ 
rium by external stresses, are common in nature and have been studied intensely over 


the last decades (see e.g. Cross & Hohenberg 

(1993 

)). Prominent examples are found in 

fluid systems (see e.g. 

Chandrasekhar ( 

1961), 

Swinney & Gollub 

(1985 

)) where the pat- 


tern formation is driven either thermally or by shear stresses. The general understanding 
of pattern forming systems has benefited from numerous experimental and theoretical 
investigations of the classical, thermally driven Rayleigh-Benard convection (RBC) in 
a layer of a simple fluid heated from below (Busse 1989 Bodenschatz, Pesch & Ahlers 
2000 [Lappa 2009). In the RBC system, the main control parameter is the Rayleigh 
number, R, a dimensionless measure of the applied temperature gradient. At a critical 
Rayleigh number, R,,, the quiescent heat-conducting basic state develops into the well- 
known periodic arrays of convection rolls characterized by a critical wavevector q^. The 
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Figure 1: Inclined convection cell of thickness d which is heated from below and cooled 
from above with temperature difference AT = Ti — T 2 > 0 for the inclination angle 
0° ^ 7 ^ 90°. Driven by gravity g the cold fluid flows downwards near the top plate and 
the hot fluid flows upwards near the bottom plate in the form of a cubic velocity profile 
(2.2). For the range 90° < 7 < 180°, the fluid layer is inverted and is heated from above. 


stability of the rolls and their evolution towards characteristic 3D patterns via sequences 
of bifurcations with increasing R has been investigated by Busse and coworkers (see e.g. 


Cross & Hohenberg (1993); Busse & Clever (1996) and references therein). The present 


paper analyzes a variant of RBC, the inclined layer convection (ILC) system, where the 
fluid layer is inclined at an angle 7 to the horizontal. Investigations of this system also 


have a long tradition (see e.g 

Vest & Arpaci 

(1969) 

Gershuni & Zhukhovitzkii 

(1969); 

Hart (1971) 

Bergholz|(1977); 

Ruth et al. 

(1980) 

Fujimura & Kelly 

(1992 

); lDaniels et al. 


p)00| )). 

In the ILC system, for 7 0° gravity g has components both perpendicular and 

parallel to the fluid layer, which leads to an important modification of the basic state 
compared to RBC. The applied temperature gradient first produces stratified fluid layers 
with continuously varying temperatures and densities. In addition, the basic state already 
contains a flow field driven by the in-plane component of g: the heavier (colder) fluid will 
flow down the incline and the lighter (warmer) fluid will flow upwards. Since the resulting 
flow field creates a velocity gradient perpendicular to the fluid layer, both buoyancy and 
shear stress driven instabilities of the basic state compete. Their relative importance 
is governed by the Prandtl number Pr, the ratio of the thermal diffusivity, k, to the 
kinematic viscosity, v, of the fluid. Furthermore, the strength of the shear stress can 
be continuously increased by increasing 7 . The orientation of the roll axes at onset of 
convection allows for directly discriminating between buoyancy and shear driving. The 
buoyancy driven rolls are aligned parallel to the incline (longitudinal rolls) while the 
shear driven rolls are aligned perpendicular to the incline (transverse rolls). It should be 
noted that the latter also bifurcate, when the fluid layer is heated from above and the 
thermal stress is therefore stabilizing. 

Our goal is not a representative parameter study of ILC with respect to R, 7 , Pr, 
which would go beyond the scope of a single paper. Our theoretical investigations have 


instead been motivated by recent ILC experiments in pressurized CO 2 (Daniels et al. 
|2000 Bodenschatz et al. 2000) with a fixed value Pr = 1.07 of the Prandtl number. In 


this work the i? — 7 parameter space has been systematically explored and a variety of 
fascinating patterns have been described. As in all ILC studies mentioned above, our 
theoretical analysis is based on the Oberbeck-Boussinesq equations (OBE). In contrast 
to the extensively studied RBC, earlier results in the literature for the ILC system are 
mostly limited to the linear regime and characterize the primary bifurcation of the con- 
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vection rolls from the basic state aX R = Rc (see 

e.g. \ 

^est & Arpaci ( 

1969 

; Gershuni 

& Zhukhovitzkii 

(1969); Hart 

(1971) 

Ruth et al. 

(1980 

; Fujimura & Kelly 

(1992)). In 

the nonlinear regime {R > Rc) Busse and Clever ( 

Clever & Busse||1977 

Busse & Clever 


1992) investigated secondary and tertiary instabilities of the convection rolls for some 


special cases. In contrast, the present work is devoted to a comprehensive theoretical 
analysis of the patterns in Daniels et al. (2000) at Pr = 1.07. 

In the present work we make use of the well-known arsenal of concepts to analyse 
pattern forming instabilities in fluid systems (see e.g. Cross & Hohenberg ( 1993| ). This 
approach deploys its full power in large aspect ratio systems (lateral extension, L, of the 
fluid layer much larger than its thickness, d) which were first realized experimentally in 
Daniels et al. ( |2000 ). From a linear instability analysis of the basic state we determine the 
critical values Rc, Qc at the onset of convection. The properties of the rolls in the weakly 
nonlinear regime, R > i?c, are then analyzed in the framework of amplitude equations, 
which yield approximate roll solutions. Using these as starting solutions allows for the 
iterative determination of the roll solutions in the nonlinear regime, where R > Rc- In a 
further step, their stability is again tested to identify the secondary instabilities of the 
roll pattern. 

We will demonstrate the agreement between experiments and theory with respect to the 
onset of convection in the 7 -i?-plane. For inclination angles 7 below a codimension 2 angle 
7 c 2 ~ 78° for Pr = 1.07 the destabilisation of the basic state is driven by longitudinal 
rolls, while transverse rolls bifurcate for 7 > 7 c 2 . Both bifurcations are always stationary 
and continuous (supercritical). The subsequent secondary destabilisation of the 2D rolls 
for 7 7 ^ 0° at increasing R is driven by oblique roll solutions, whose axes are not along 
the longitudinal or transverse roll directions. As a result, spatially periodic 3D patterns 
are often observed. These are characterized by the nonlinear interaction of three roll 
modes with wave vectors ( 51,92 and 53), that fulfil a wavevector resonance condition 
5i + 52 + 53 = 0. 

As common in other large aspect ratio convection experiments, one also finds imper¬ 
fectly periodic, weakly turbulent patterns. For instance, the 3D motifs mentioned above 
appear locally superimposed on the original the 2D roll pattern, where they burst and 
vanish repeatedly in time (Daniels et a/. [2000 Daniels & Bodenschatz 2002 Daniels et al. 


[2003 ). To test whether such dynamic states are caused by experimental imperfections 
(e.g. lateral boundaries, spatial variations of the cell thickness, inclination of the cell in 
two directions), we have performed the first comprehensive numerical simulations of the 
OBE in ILC for large aspect ratio convection cells. While a conclusive assessment of the 
underlying mechanism producing bursts remains elusive, the weakly turbulent dynamics 
of the pattern have been well reproduced in our simulations. Since the shear stresses play 
an important role in our system, there might be an analogy to the instabilities of the 
laminar state in purely shear driven fluid systems like Couette or Poisseuille flow, which 
also often appear in the form of localized events (for recent examples see e.g. (Lemoult 


et al. 2014 Tuckerman et al. 2014)). We hope that investigations such as this one will 


reveal deeper commonalities between ILC and such shear driven patterns in the future. 

The paper is structured as follows: a brief summary of the governing OBE for the 
ILC system is given in We then discuss the onset of convection in terms of 5 c, Rc 
for Pr = 1.07 and the resulting periodic roll pattern in The results of the stability 
analysis of the rolls the nonlinear regime and the resulting phase diagram are presented 
in In (J^we show direct simulations of the OBE for different R and 7 , which compare 
well with the experiments. A short summary of this work together with perspectives for 
future work can be found in In three detailed appendices we present first the detailed 
OBE equations for the ILC system and discuss the numerical method to characterize 
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the roll solutions above onset together with their secondary instabilities (Appendix [A|) . 
Then we address briefly our approach to solve the OBE in general using direct numerical 
simulations (Appendix]^. Finally we return to the linear stability analysis of the basic 
state to give additional informations regarding the properties of qc, Rc (Appendix [C|) . 


2. Oberbeck-Boussines equations for ILC 

As shown in figure we consider convection in a fluid layer of thickness d, which 
is inclined at an angle 7 (0° < 7 < 180°) with respect to the horizontal. Constant 
temperatures T 2 (Ti) with difference AT = Ti — T 2 >0 are prescribed at the boundaries 
(z = ±d/2) of the layer. Both cases, heating from below (0° < 7 ^ 90°) and heating 
from above (90° < 7 < 180°), will be considered in this paper. 

The resulting ILC system is described by the standard Oberbeck-Boussinesq equations 
(OBE) for incompressible fluids. As usual, the OBE are non-dimensionalized using d as 
the length scale and the vertical diffusion time = d ?/k as the time scale. The velocity 
u is measured in units of d/t^ and the temperature T in units of Tg = vn/agd^ with a 
the thermal expansion coefficient. Using a Cartesian coordinate system aligned with the 
layer (see figure [^, the OBE read as follows: 

[d/dt+{u-V)]T = V‘^T PRz-u, (2.1o) 

Pr-^ [d/dt + (m • V)] M = -^T-Vp, (2.16) 


where V • m = 0 due to incompressibility and g = —g (cos 72 + sin 7 £) describes the 
effect of gravity with the gravitational constant g. All terms which can be expressed as 
gradients are included in the pressure term Vp. Equations (2.1 1 are characterized by the 
angle of inclination 7 along with two nondimensional parameters, the Prandtl number 
Pr = VIK, and the Rayleigh number R = AT/Tg. 

In line with previous theoretical investigations of ILC in the literature (see in particular 
(Clever & Busse 1977 Busse & Clever 19921), we idealize the system to be quasi-infinite in 
the x — y plane. This is considered to be the appropriate description for large aspect-ratio 
systems. Equations (2.1) then admit primary (basic) solutions (denoted with subscript 
0) of a linear temperature profile To(z) and cubic shear velocity profile Uo{z): 


\Ti+T2 1 


\ 2 1 ] 

2AT 

, Uo[z) =xsm^R- 

r ■ 4J 


X sin-fRU^iz). 


( 2 . 2 ) 


It is convenient to introduce modifications 9 and v of the basic state and describe the 
secondary convective state as: 


T{x,z,t) =To{z) + 9{x,z,t), u{x,z,t) = Uo + v{x,z), X = {x,y), (2.3) 

which fulfill the boundary conditions 9{z = ±1/2) = v{z = ±1/2) = 0. Furthermore, the 
solenoidal velocity field v is mapped by the well-known poloidal-toroidal decomposition 
to two scalar velocity functions /, $(a;, z, t) and a correction U (z, t) of Uo{z); for details, 
see Appendix The resulting coupled set of equations for 9, /, 4), U are analyzed in 
the following sections using standard Galerkin methods and direct numerical simulations 
(DNS). 


3. Finite-amplitude roll solutions 

Spatially periodic convection roll solutions of the OBE (2.1) with wavevector q exist for 
Rayleigh numbers R > Rc, where the homogeneous basic state (2.2) is unstable against 
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infinitesimal perturbations which depend on x, y. The onset of convection in ILC system 
at the critical Rayleigh number associated with the critical wavevector has been 


Gershuni & Zhukhovitzkii 

(1969) 

Birikh et al. 

(1972 

); |Hart 

(1971 


useful overview can be found in Chen & Pearlstein (19891 and references therein. Some 
additional general information is given in Appendix [Cj 

Since the ILC system is anisotropic, we have to consider the linear stability of the basic 
state against arbitrarily-oriented convection rolls with wavenumbers q = (/(cos' 0 , sin ■)/)). 


For that purpose, we have to analyse (2.1) linearized about the basic state (2.2). For 


details of the standard numerical method, see Appendix A.l For 7 = 0 (horizontal layer 
with Uq = 0) the system is isotropic and we have the standard Rayleigh-Benard convec¬ 
tion (RBC) where \qc\ = Pco = 3.1163 and Rc = Rco = 1707.762 (see e.g. Dominguez- 


Lerma et al. (1984)) which depend on neither ip nor Pr. This is distinct from finite 7 , 


linear equations (see A4&). 


since Uo{z) defined in (2.2) yields a contribution proportional to cos"0siny/Rr in the 


In the following, we concentrate on the special case Pr = 1.07, where the bifurcation of 
the basic state is always stationary; other Pr are briefly discussed In Appendix |C.2| Figure 
[^displays the rescaled critical Rayleigh number Rc/Rco and the critical wavenumber 
as function of the inclination angle 7 and different ip. In general, only two particular 
q—orientations turn out to be relevant (see e.g. Appendix |C.2[ ). The convection solutions 
at onset are either buoyancy driven longitudinal rolls with their axes along the incline, 
i.e. qc = PcV, Ip = 90° or shear driven transverse rolls with their axis perpendicular to 
the incline where qc = QcX, ip = 0°. 

Longitudinal rolls {ip = 90°) exist only in the range 0 < 7 < 90° (heating from 
below, see figure . Their critical wavenumber is given by (/(, = qco for all 7 and the 


critical Rayleigh number i?(,( 7 ) fulfils the relation i?(.( 7 ) cosy = i?cO (see Appendix C.2), 


implying that i?(, diverges in the limit 7 —>■ 90°. In contrast, a bifurcation to transverse 
rolls exists in the whole interval 0 < 7 < 180°. The critical Rayleigh number A* ( 7 ) rises 
continuously as function of 7 and diverges at 7 = 180° (stable horizontal fluid layer, 
heated from above). In figure]^ the critical data have been shown only for 7 up to 120°, 
where Rp ~ 10^ involves large thermal gradients. Thus, the use of the OBE becomes 
questionable for 7 > 120°, since non-Boussinesq effects due to temperature variation of 
the various material parameters should be taken into account. 

Inspection of figure reveals the existence of a codimension-2 bifurcation point 7 c 2 = 
77.746° where = i?* = 8046.420, such that for 7 < 7 c 2 longitudinal rolls bifurcate at 


onset { RI < Rl) while for ' 

i > 7 c 2 the transverse ones prevail, i 

in Gershuni & Zhukhovitzkii 

-TT"?——- 1 tt: -' 

(1969) and detailed in Appendix C.2 


tions of the critical values Rcil) and (/‘(y) of the transverse rolls. In this paper, we will 
often use the reduced main control parameter e defined as: 


e={R-Rc{j,Pr))/Rc{j,Pr) 


(3.1) 


as a measure for the relative distance from threshold Rcij) at e = 0, instead of R. 

The standard computational methods to construct finite-amplitude roll solutions with 
wavevector qc for R > Rc, where exponential growth of the linear modes is balanced 
by the nonlinear terms in the OBE, are sketched in Appendix |A.2| The amplitudes of 
the roll solutions grow continuously like -^/e (see the discussion after (A 8 )). Thus, the 
primary bifurcation to rolls is continuous (forward). 
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Figure 2: Normalized critical Rayleigh number Rc/Rco (left panel) and critical wave 
number \qc\ (right panel) as function of the inclination angle 7 for Pr = 1.07 for different 
roll orientations ip. Dotted line is for ip = 90° (longitudinal rolls), dashed line for tp = 60°, 
dashed-dotted line for tp = 45° and solid line for '0 = 0° (transverse rolls). 


4. Secondary instabilities of roll solutions for Pr= 1.07 

In this section, we discuss the secondary destabilisation mechanism of rolls with wavevec- 
tor Qc, in a ILC system with Pr = 1.07 and inclination angle 7 , that become unstable 
when e = emstil)- Based on the methods described in Appendix A.2 the stability dia¬ 
gram presented in figure has been determined in the 7 — e plane.The solid lines mark 
the locations of the various secondary instabilities of the finite amplitude roll solutions 
with q — qd'j) at e = einst[l)- Thus to the left of this line for 7 « 15° and below this line 
for larger 7 stable roll solutions exist. For details of the dependence of these secondary 
stability lines upon cut-off parameters in the Fourier-Galerkin expansions see Appendix 

According to figure the type of secondary roll instabilities depends strongly on the 
inclination angle 7 . Our main interest in this paper are the various ZD patterns which 
develop for e > Cmst- As a first impression, we include in figure excerpts of patterns 
observed in experiments ( Daniels et aLpOOO l. Here, we aim to reproduce and interpret 
such 3D patterns based on direct numerical simulations (DNS) of the OBE. For details 
of time stepping scheme we refer to the algorithm presented in Appendix [B| 

We first consider small e < Cinsdl)- Starting from random initial conditions, modes 
with wave vector qd^i) prevail leading to perfect roll patterns, as shown in figure 
The DNS are performed on a square with side lengths = Ly = 12Ac with Ac = 
2'K/qcdi) where we obtain, longitudinal rolls at 7 = 10 ° and transverse rolls at 7 = 
85°. Here, and in the rest of this paper, we show snapshots of the vertical {z) average 
{9{xd)) of the temperature field 9{x,z,t) (see (A5)). Throughout this paper the height 
of the convection cell increases from left to right, i.e. with increasing x with respect 
to the coordinate system attached to the cell (see figure [^. Such pictures, which we 
refer to as the temperature plots, are typically used in the literature to compare with 
experimental convection patterns, which are visualized via shadowgraphy (for examples, 
see Bodenschatz et al. (2000)). In gas convection experiments of the type considered 
in this paper, dark and bright regions in figure indicate positive (hot) and negative 
(cold) variations in the temperature field around the basic linear temperature profile 
( de Bruyn et qZ.[ 19M ). We use a 8 -bit grayscale to visualize {9{x)); whose range increases 
monotonically as a function of e. 
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Figure 3: Phase diagram of the convective roll patterns in the 7 ,e plane for Pr = 1.07. 
Solid lines (coloured online) indicate the secondary instabilities of the primary roll pat¬ 
terns at e = Cinst, he. at i? = {l-\-einst)Rci'y)- The codimension 2 point with jc 2 = 77.746° 
is marked by a star. For 7 < 7 ^ 2 , we obtain longitudinal rolls (LR), and for 7 > 7 c 2 , 
transverse rolls (TR). Increasing 7 in the interval 0 ^ 7 < 7 c 2 , we obtain the thresholds 
for skewed varicose instability (SV) shown in purple, longitudinal subharmonic oscilla¬ 
tions (LSO) shown in light blue and wavy rolls (WR) shown in dark red as instabilities 
of the LR. In the range 7 > jc 2 the thresholds for instability of the transverse rolls 
through the knot instability (KN) is shown in black and for the transverse oscillations 
(TO) in brown. The graph is decorated with excerpts from corresponding experimental 
pictures (Daniels et al. 2000| . In addition crawling rolls (CR), transverse (TB) and the 
longitudinal bursts (LB) are shown, which cannot be directly associated to the secondary 
instabilities. Black open circles indicate locations in the 7 — e plane, where experiments 
and numerical simulations are compared in section 


In the following, we discuss the secondary instabilities of the primary convection rolls 
in detail. We start with inclination angles 7 below 7 c 2 in §4.I[ before we concentrate on 
the vicinity of 7 c 2 and finally on larger 7 . 

The various secondary instabilities are visualized by direct simulations of the underly¬ 
ing OBE (see Appendix |B.I[ ) . In general we use a minimal rectangular integration domain 
in the x — y plane which is consistent with and the wavevectors of the dominant desta¬ 
bilizing modes. For visualisation, the domain is periodically extended to a larger domain 
with Lx = I2Ac = Ly. 


4.1. Secondary roll instabilities below 7 c 2 

In this section, we will characterize in detail the secondary instabilities of longitudinal 
rolls in figure 


4.I.I. Skewed varicose instability (SV) 

For small inclinations (7 ^ 5°), we recover the well known skewed varicose (SV) 
instability for planar RBC with 7 = 0 ( Busse fc Clever||l979 1. This is a stationary long- 
wavelength instability where the original longitudinal rolls are slowly modulated along 
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Figure 4: Temperature plots {9{x)), for Pr = 1.07 from direct numerical simulations of 
the OBE (2.1): (a) buoyancy dominated longitudinal rolls oriented along the incline from 
right to left at Pr = 1.07, e = 0.01, 7 = 10° with = (0,3.1163) (left panel); (b) shear 
dominated transverse rolls oriented perpendicular to the incline at e = 0.02, 7 = 85° 
with Qc = (2.82,0) (right panel). In the plane the figures are oriented along the x-axis 
of local coordinate system in figure(reproduced at the right). All temperature plots in 
the rest of the paper are shown using the same convention. 




Figure 5: Subharmonic oscillatory instability of longitudinal rolls (LSO) for Pr = 1.07, 
7 = 17° and e = 1.3: temperature plot (left panel) together with the wavevectors qi = 
Qc, < 72,93 (wavevector resonance <71 = 92 + 93 ) of the leading Fourier modes (right panel). 


their axes but also with respect to their distance. The SV instability will not be further 
discussed in this paper. 

4.1.2. Logitudinal subharmonic oscillations (LSO) 

In the range 5° < 7 < 21°, the longitudinal roll pattern with q = q^ = (0,3.1163) 
becomes linearly unstable to oscillatory subharmonic perturbations with wavevectors 
92,3 = (+ 9 a;, 9 c/ 2 ) and a finite circular frequency ujinst- For the representative case 7 = 
17°, where primary rolls get unstable at Cinst = 1-044, we have = 1.279 and = 
10.21. In figure we show an excerpt from our simulation at e = 1.3. performed on 
the minimal rectangle with sidelengths = ^'njqx and L^ = 2 Ac/ 9 c with Ac = 2t: jqc. 
The periodically extended picture is six times larger. The pattern is characterized by 
periodic modulations of the longitudinal rolls, which are in phase on every second roll 
reflecting the subharmonic nature of the instability. The resulting LSO pattern obeys the 
wavevector resonance 91 — 92 — 93 = 0 as seen in the Fourier spectrum. 

The time evolution of the pattern in the x — yplane is: 

f{x, y) = Acos( 9 c 2 /) + B(t) cos^q^x + qcy/2) + C{t) siniq^x - qcy/2) (4.1) 
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Figure 6 : The coefficients A (horizontal upper dotted line), B(t) (dash-dotted) and C(t) 
(solid line) in (4.1) as function of time in units of T^; B(t), C(t) are multiplied by a factor 
of four for better visualisation. 


The coefficient A is independent of time while B(t) and C(t) are periodic with circular 
frequency uj « Wmst- The time evolution of these coefficients (given in units of Tg = 
0(lmK) defined in are shown in figure]^ 

The LSO instability has been first described by Busse & Clever ( |2000 ) at the lower 
Prandtl numbers Pr = 0.7, by solving ( |A7t in the minimal integration domain. In 
addition for certain parameter combinations of Qx , Qy an intermittent appearance of bursts 
has been described, which we have not reproduced for Pr = 1.07. 

4.1.3. Wavy Roll (WR) instabilities 

The next instability type, characterized by the appearance of longitudinal rolls un¬ 


dulating like snakes along their axes, has been first described in Clever & Busse (19771 


where the notion wavy-instability has been coined. The resulting bifurcation to wavy 
rolls (WR) is observed in a fairly large y-interval between 21° < 7 ^ 7 c 2 , very close to 
onset of convection ( Cmst = 0(0.01)). In the framework of the Galerkin stability analysis 
in Appendix |A.2[ this instability is characterized by long-wavelength destabilizing modes 

with wavevectors Qmax = i±qx,qc) with (|ga;| < [gd)- _ 

The WR have been discussed in detail in [Daniels et al. (20081, to which we refer 
for details. Here one finds representative experimental and theoretical pictures (see also 
figure [T^ in )|^below) as well as a Galerkin stability analysis of the 3D wavy roll patterns. 
In fact, for e > Cmst stable WR with finite qx exist. They are spatially periodic in the 
plane, characterized by the wavevector resonance 92 + <73 = 2qc with < 72,3 = (±<7a:) <?c)- In 
the X, j/—plane the temperature pattern is described as 

f{x,y) = Acos{qcy) B sm{qcy) sm{qxx). (4.2) 

For e = 0.1 and qx = 1.28, we find for instance A = 48.1 and B = 12.8 in units of Tg. 


4.2. Secondary roll instabilities above 7 c 2 

In this section we discuss the secondary instabilities of the transverse rolls bifurcating 
for inclinations 7 > 7 c 2 - 
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Figure 7: Knot instability of transverse rolls at Pr = 1.07, 7 = 81.9° and e = 0.055: 
temperature plot (left panel); wavevectors of the leading Fourier modes gi = Qc, 92 = 
( 0 , 9 co) and qa = qi + q2 (right panel). 


4.2.1. Knot (KN) instability 

Just above the codimension 2 point, the shear dominated transverse rolls with qi = 
{qc, 0 ) and ~ 2.82 (see figure]^ are destabilized by the longitudinal rolls with wavevec- 
tor 92 = (0; %) and Qy = Qcq = 3.1163. The steeply rising stability line starts at 7 = 7 c 2 
at e = 0. In view of the logarithmic e-scale, the curve is only shown for e > 0.01. In 
the weakly nonlinear regime, the oblique mode with wave vector 93 = ( 9 ^, 9 ^) is already 
important and the wavevector resonance 93 = 9 i -b 92 is established. The pattern is well 
described by: 


f{x, y) = A cos{qcx) -b B sui{qyy) + C sin( 9 ca:) cos{qyy) 


(4.3) 


As a function of e, the amplitudes B,C increase continuously above e = Cinst- 

In figurej^we address an representative example for 7 = 81.9° with Emst = 0.053, where 
originally Lx = Ac and Ly = 2Trjqco. At e = 0.055, we find A = 17.9, B = 1.95, C = 2.16 
in ( |4.3| ) . The resulting stationary pattern (see fig ure [7| left panel) has some similarity 
to the knot patterns described in iBusse & Clever! (|1979|) for the isotropic RBC system. 


However, contrary to the patterns in Busse & Clever (19791, we observe a wavevector 


resonance triggered by the oblique mode 93 in the ILC system. 

For completeness, it should be mentioned that the knot instability in ILC has been 
previously investigated in Fujimura & Kelly (1992) for 7 < 90° in the framework of two 


coupled amplitude equations restriced to the amplitudes A,B in (4.3). 


4.2.2. Transverse Oscillatory rolls (TO) 

For 7 > 83.2° the destabilization of the transverse rolls starts to be governed by the 
transverse oscillatory rolls (TO) along an almost horizontal transition line as function 
of 7 . Transverse oscillatory rolls (TO) are characterized by destabilizing modes with a 
Floquet vector s of relatively small but finite modulus |s| ~ 9 c /6 and by an oscillatory 
time dependence of period about 3.5t„. In an analogy to the stationary SV instability 
of longitudinal rolls (see (4.1.1), the rolls are expected to become slowly modulated 
along their axis and also with respect to other in-plane direction. Without a definitive 
resonance condition among the dominant destabilizing modes, we do not expect periodic 
3D patterns of the kind discussed in the previous subsections. Thus, we have to solve 
the OBE on a larger domain in the x — y plane, including modes with wave vectors 
|(9x,9y)| ^ 9c- An excerpt of representative DNS pattern for 7 = 84.9° is shown in 
figure 1^ (left panel) at e = 0.07. This pattern results from a secondary instability of 
transverse rolls (i?c = 8282.64 and q^ = 2.8023) at Cinst = 0.063 with Winst = 1-809. 
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Figure 8 : Transverse oscillatory rolls at Pr = 1.07, 7 = 84.9° and e = 0.07: temperature 
plot (left panel); wave vectors of the dominant Fourier modes (right panel). In this case, 
we use a square grid with Qy = Qc = 3.117. 


One observes the appearance of localized patches with reduced amplitudes on top of 
the original slowly modulated transverse rolls. Apparently, the pattern arises through a 
complex beating phenomena, due to a superposition of oscillating modes with slightly 
different wave vectors. The most prominent ones ( 92 , Qs) together with qi = are shown 
in figure(right panel). 

Figurel^shows the complicated evolution of the pattern during one cycle {2tt/ uinst)- 
The localized patches of reduced amplitudes in figure [^i), evolve into slanted lines of 
reduces amplitudes in (ii). Further in the cycle, the undulation amplitude of the rolls 
first increases and decreases then before arriving again at the initial pattern in figure 


l^iii-vi). 

According to our phase diagram in figure the instability of the transverse rolls 
towards the TO pattern remains relevant for larger 7 and also governs the secondary in¬ 
stability in the case of heating from above. A representative example of a time sequence is 
shown in figurefor 7 = 100°, where the TO instability is characterized by Cinst = 0.06 
and ujinst = 1.48. The graph of the most relevant destabilizing modes looks practically 
identical to the one in figure|^and is thus not shown. The interaction of the modes leads, 
however, to a much simpler time evolution, as compared to figure]^ for 7 = 84.9°. It is 
possible that in the latter case, the destabilizing modes triggering the knot instability for 
slightly smaller 7 in hgurej^come into play as well. 

As will be discussed in the regions of suppressed amplitudes become elongated 
and are no longer periodically arranged in the plane for increasing aspect ratio {Lx,Ly). 
As shown i n figure |18| bel ow, they compare well with the corresponding experimental 
patterns in (Daniels et a7|[2000), called switching diamond panes (SDP) there. 


4.2.3. Vertical convection 

The case of a vertical convection cell (7 = 90°) is of special interest since the pattern 
formation is exclusively driven by the shear stress. So this system has motivated many 
previous investigations, mainly in the linear regime and often with Pr > 12.45 where an 
oscillatory bifurcation to transverse rolls takes place. For Pr = 1.07, however, we observe 
stationary transverse rolls at onset, and their stability analysis yields always a secondary 
bifurcation to the TO pattern for all 7 near 90° (see figure and our discussion in the 
previous subsection). 

This result is noteworthy, since in the previous literature ( [Clever fc Busse 19951 for 
Pr = 0.71 (air) a stationary secondary instability of the transverse rolls driven by the 
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Figure 9: A time sequence of TO patterns shown over one time period for the case of 
heating from below. Consecutive panels are separated by QAAt^. System parameters are 
the same as figurej^with Pr = 1.07, 7 = 84.9° and e = 0.07. 



Figure 10: A time sequence of TO patterns shown over one time period for the case of 
heating from above. Consecutive panels are separated by QAt^. System parameters are 
Pr = 1.07, 7 = 100° and e = 0.08. 


effective subharmonic roll modes with wavevectors < 72,3 = (9c/2,±p) was predicted. This 
finding has been confirmed by our own calculations. The instability of the transverse rolls 
with Qc = 2.8123, i?c = 5701.2625 takes place at einst = 0.0599 with p = 1.5898. These 
numbers are consistent with those used for direct simulations of the OBE in IClever fcl 
|Busse| ( [19^ (i?c = 5726.9, (7c = 2 69,p = 1.7 and e > 0.11). 

In close analogy to the LSO (14.1.2), the resonance conditions <7c = <72 + <73 holds. 
To confirm the results of stability analysis we have performed DNS of the OBE on the 
smallest periodicity domain in the plane compatible with instability data above, i.e. 
with Lx = 2\c,Ly = Xc{qc/p)- The resulting stationary temperature plot (periodically 
extended) for e = 0.064 > einst is shown in figure 11 for = Ly = 12Ac. Keeping the 


dominant Fourier modes, the pattern is represented as: 


f{x, y) = A cos{qcx) - B sm{qcx/2) cos{py) (4.4) 

with A = 131.23, i? = 60.76 in units of T^; the subharmonic instability of the transverse 
rolls is obvious from the argument of sine in the second term. 

A closer look at the Pr-dependence of the secondary bifuraction of the transverse rolls 
in this regime shows that for Pr > 0.9 the secondary SFIV bifurcation of the transverse 
rolls is replaced by the TO bifurcation discussed in §4.2.2[ This is consistent with our 
stability analysis and experimental observations at Pr = 1.07. 
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Figure 11: Subharmonic varicose instability of transverse rolls with Pr = 0.71, e = 0.065 
and 7 = 90°: temperature plot (left panel); wavevectors of the leading Fourier amplitudes 
(right panel). Given that Qx = 2.8123 and Qy = 1.59 and choosing the number of rolls 
ul = 12, we have Aqx = qxjnh = 0.2344 and = qy/uL = 0.1358. 


5. Comparison with experimental results 

In the previous section, we have discussed the various characteristic secondary insta¬ 
bilities of the ILC roll patterns with increasing inclination angle 7 . A number of basic 
destabilization mechanisms have been identified by considering simulations in small pe¬ 
riodicity domains in the plane of linear dimension L, where L = 0{2d) with d the cell 
thickness. Our goal in this section is a comparison with the pressurized CO 2 experiments 
( Daniels et aZ.|[2000 Daniels|[2002[ ) at Pr = 1.07. In these experiments a convection cell 
with very small layer thickness d = (710 ± 7)p, m could be realized together with quite 
large lateral dimensions [(42 x 21)d^]. Thus, the dimensions of the convection cell are 
such that the experiments are expected to be well described in simulations by periodic 
boundary conditions in the plane. Furthermore, such dimensions give a vertical diffusion 
time ty = d?/K = 3.0 s, which sets the time scale, is small enough for typical experiments. 
In the rest of this section, we will present the generic features of the patterns when all 
transients have died out. To compare with theory, we have performed numerical simula¬ 
tions of the basic equations ( 2 . 1 ) on a large horizontal domain with lateral dimension up 
to L = 20Ac = O(40d). For the numerical method, see Appendix 

Excerpts of the experimental shadowgraph pictures to be discussed in this section are 
already shown in figure We will follow the same sequence of parameter combinations 
used in the preceding section, where 7 was systematically increased. As discussed in 
we show in this section the vertical temperature average {9{x)) side by side with the 
corresponding experiments. A quantitative agreement between theory and experiment 
is not to be expected, as along with the complicated optics involved in shadowgraphy 
( Trainoff fc Canell||2002 1 , the experimental pictures are typically digitally remastered to 
enhance their contrast. 


5.1. Convection pattern for 7 < 7 c 2 

We start with the subharmonic oscillatory patterns (LSO), in figure which bifurcate 
from longitudinal rolls. As shown in figure [T^ experiments and theory match very well. 
We observe patches of subharmonic oscillations, which we have previously discussed in 
§4.1.2| on the basis of a stability analysis of the longitudinal rolls and numerical simula¬ 
tions in a minimal domain (containing only one roll pair). In large ILC systems, typically 
such motifs appear only in localized patches that compete with moderately distorted 
rolls. Such patches expand and shrink in time and their centers move erratically over 
the plane. The subharmonic oscillations within the localized patches show an internal 
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Figure 12: Snapshots of an longitudinal subharmonic oscillations (LSO) for Pr = 1.07, 
7 = 17° and e = 1.5 from experiments (left panel) and from our numerical simulation 
(right panel). 




Figure 13: Representative wavy roll patterns (WR) from experiments (left panel) in 
comparison with simulations (right panel) for Pr = 1.07, 7 = 30° and e = 0.08. 


dynamics with a time scale of 1 to 3 cycles per which is of the same order as the 
period jiriinst of the Hopf bifurcation in ( 4.1.2| 

With increasing 7 , the longitudinal rolls become unstable against undulations even 
for very small e; i.e. the instability line in figure bends dramatically down. Typical 
experimental and theoretical pictures of undulated (wavy) rolls (WR) shown in figure [T^ 
again match very well with each other. However, instead of the stationary undulations 
as predicted in §4.1.3[ the patterns are characterized by patches of uniform undulated 
rolls, separated by grain boundaries ( Daniels et a/.|[2000 Daniels fc Bodenschatz|[2002 ). 
In addition, the rolls are scattered with point defects that move at right angles to the 
rolls. The wavy patterns have been discussed in detail in Daniels et al. (2008), where also 
a weakly chaotic dynamics of the amplitudes A,B in (4.2) is analyzed in detail. 

With increasing e, the undulations become more and more disordered and the rolls get 
disrupted. A transition is observed to the dynamic state of the so called crawling rolls 
(CR) ||Daniels et al. (2000) as shown in figure 14 This state is also obtained in our nu¬ 
merical simulations, which indicates that it is not caused by experimental imperfections. 


5.2. Convection close to codimension 2 point 
The vicinity of the codimension 2 point 7 c 2 is of particular interest. According to figure 
1 ^ the wavy roll instability governs the secondary instability of the longitudinal rolls 
up to 7 = 7 c 2 - In contrast, for 7 > 7 c 2 the primary transverse rolls are predicted to 
become unstable against cross rolls leading to the knot patterns (KN). This instability 
mechanism is confirmed by the pictures shown in figures |15[ It is remarkable that the 
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Figure 14: Crawling rolls (CR) from experiments (left panel) and from simulations (right 
panel) at the same parameters as in figure 


13 


except at a larger e = 0.4. 



Figure 15: Knot pattern (KN) at 7 = 80° and e = 0.05 in experiments (left panel) and 
in simulations (right panel). 



Figure 16: Localized transverse bursts (TB) for the parameters Pr = 1.07, 7 = 77° and 
e = 0.04 in experiments (left panel) and in simulations (right panel). 



DNS on a large domain in the plane {L^ = Ly = 20Ac) starting from random initial 
conditions has led to perfect knot patterns. They are indistinguishable from those shown 
in figure generated on a small domain with Lx,Ly « Ac. The teeth-like structure on 
the transverse rolls caused by a resonant interaction of the three roll modes (see figure]^ 
is born out in the experimental picture. However, the transverse rolls are slightly oblique 
here and undulated. 

We now discuss two types of patterns which do not allow for a direct interpretation 
by secondary instabilities of the basic rolls. First, we show in figure the transverse 
bursts (TB) for 7 = 77° < 7 c 2 and at e = 0.04 slightly above the secondary wavy 


























































Figure 17: Localized longitudinal bursts (LB) for the parameters Pr = 1.07, 7 = 79° and 
e = 0.1 in experiments (left panel) and in simulations (right panel) 



Figure 18: Switching diamond panes (SDP) for the parameters Pr = 1.07, 7 = 100° and 
e = 0.1 in experiments (left panel) and in simulations (right panel). 


bifurcation of the longitudinal rolls. Both in experiments and simulations, we observe a 
background of slightly undulated rolls with some amplitude modulations. Intermittently, 
localized transverse structures (bursts) appear, which contract, vanish and reappear at 


other places. The longitudinal bursts in experiments have been analyzed in Daniels et al. 


(2003), to which we refer for more details. 


In contrast, for 7 > 7 c 2 and intermediate e, longitudinal bursts (LB) are observed; rep¬ 
resentative examples are shown in figure [P^ They are characterized by localized loops of 
longitudinal rolls superimposed on transverse rolls. The experimental picture shows more 
of the bimodal knot pattern in the background than do the simulations. In the vicinity 
of 7 c 2 5 we do not expect simulations to reproduce all details of the the experiments at 
the same parameters, as the system is very sensitive against small changes of 7 and e. 
The material parameters in the experiments certainly have some inaccuracies. In addi¬ 
tion, non-Boussinesq effects presumably lead to a slow drift of the experimental pattern. 
Although the two burst phenomena are clearly reflected in our simulations, additional 
efforts are necessary in the future to understand their underlying mechanism. 

5.3. Shear stress dominated instabilities 

Finally, we briefly address the heating-from-above case which is described in for which 
the inclination angle is 7 > 90°. As discussed before, the destabilization of the basic state 
is due to the shear stress of the cubic flow profile Uq (2.21. In ( |4.2.2 we have described 
the destabilisation of the primary transverse rolls to switching diamond pane patterns 
shown in Daniels & Bodenschatz (2002). According to figure]^ the transition line to 
the TO is almost horizontal and begins at 7 slightly above 7 c 2 - In figure [T^ we show a 
































Figure 19: Patterns formed during a typical evolution of switching oscillatory rolls shown 
in order, for the case of heating from above as in figure |18| Consecutive panels are 
separated by 5ti,. System parameters are Pr = 1.07, 7 = 100° and e = 0.1. 



Figure 20: Chaotic switching rolls at 7 = 100°,e = 0.19 [experiments (left panel), simu¬ 
lations (right panel)]. 

representative example for 7 = 100 °, where experiment and simulations agree very well. 
The time evolution observed in the corresponding time sequence, presented in figure [T^ 
reflects the frequency uJmst = 1-48 given in §4.2.2| and documented in figure[^ Increasing 
e further causes the patches with enhanced amplitude to get smaller and to move more 
erratically as shown in figure | 20 | 


6. Conclusions 


The recent experimental study of ILC for Pr = 1.07 by Daniels et al. (2000) has 


opened a new path to a much better understanding of this dynamically rich system. 
For a recent overview, see chapter 7 in Lappa ( 2009[ ). In contrast to previous work, the 
convection instabilities of the basic state, sketched in figure have been systematically 
explored as function of the inclination angle 7 and the Rayleigh number R. Furthermore, 
the resulting patterns are visualized directly and the large variety of new pattern types 
are shown in the phase diagram in figure]^ 

For the theoretical analysis in this paper, the large lateral extent of the convection cell 
in the experiment has been of particular importance. For large aspect ratio systems the 
influence of the lateral boundaries of the cell is certainly highly suppressed, making it ap¬ 
propriate to use periodic boundary conditions. A convincing agreement with experiments 


by Daniels et al. (2000) has been obtained. The analysis of the primary bifurcation of 


rolls at the onset of convection and their secondary bifurcations have revealed the compli¬ 
cated interaction of buoyancy and shear driven destabilization mechanisms. Of particular 
importance is the spatially resonant interaction between three roll systems with different 
orientations in the plane (wavevector resonance, as detailed in Q. 

A look at the experimental pictures and the pattern dynamics (see 


Daniels et al. 
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(2000)) shows that they are not completely described by perfect periodic patterns either 
in one or two dimensions in the fluid layer plane. One finds cases where a kind of clear 
periodicity is expressed only in parts of the cell (figures . One also observes defect 

lines; in addition the patterns change in time. Further increasing R in these cases often 
leads to turbulent patterns (see figures [20| ) . In addition, there are other cases where 
localized patches of a different structure than the underlying, regular background pat¬ 
terns appear intermittently. Two examples are the transverse bursts in figure and the 
longitudinal bursts in figure The ability to reproduce such weakly turbulent patterns 
in direct numerical simulations of the OBE validates their generic character, independent 
of experimental conditions. 

To unravel the basic underlying mechanism is a difficult task; here we have been 
unable to understand these states in terms of instabilities of the underlying roll patterns. 
Phenomena where complex patterns that cannot be explained in terms of the instabilities 
of the underlying simpler patterns, have been previously described in standard RBC. A 
prominent example is the spiral defect chaos (Bodenschatz et al. 2000), which is often 
observed for medium Pr and for Rayleigh numbers R slightly larger than Rj., where rolls 
are linearly stable. 

Further effort is thus needed to analyze and to quantify the dynamics of the turbulent 


events in detail as has been done for the bursts in Daniels & Bodenschatz (2002); Daniels 


et al. (2003) or for the wavy patterns in Daniels et al. (2008). Another issue is the weakly 
turbulent convection states described by Busse and coworkers in ILC, appearing even 


for very small systems containing only one roll pair ( Busse fc Clever||l992 2000). Their 
relation to the weakly turbulent events, which here cover considerably larger domains in 
the plane, requires further investigation. 

In this paper, we have restricted ourselves to the special case of Pr = 1.07. As part 
of future work, it is planned to apply our methods in particular to fluids with large 
Pr > 12.47, where the primary roll bifurcation is oscillatory. 


The authors are highly indebted to Prof. F. Busse for his very useful comments and 
fruitful discussions on the subject of this paper. 


Appendix A. Governing eqnations and stability of rolls 

In section the poloidal-toroidal decomposition of the solenoidal velocity field v (2.3) 
is written as follows: 


v{x,y,z,t) = V X (V X fz) -I- V x $£ -I- U{z,t) = xf + + U{z,t). 


(Al) 


The explicit equations for 0,/, 4) are obtained by inserting (2.3) into (13 followed by 
requiring the divergence of velocity to vanish in (2.16). The evolution equation for the 


secondary meanflow flow U{z,t) is obtained by averaging the velocity equation (2.16) 
over the x — y plane, leading to: 


1 dU{z,t) 1 d{v,v) d'^U 


+ -bsiny^- {dx,dy,0){PxX + Pyy), 


(A 2) 


Pr dt Pr dz dz^ 

where the overbar indicates a horizontal average. 

Except for minor changes, the resulting equations for 9, /, $ can already be found 
Busse & Clever (1992); Daniels et al. (2008|). The differences arise firstly from the 


definition of the Rayleigh number R = AT cosy/T^, whose explicit dependence on cosy, 
is not convenient for the description of vertical convection cells when y = 90°. By the 
transformation 0—^0/cosy and interchanging x and y we arrive at our formulation. 
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Secondly, our equations for the meanflow U (A 2) contain the additional pressure terms 
Px{t),Py{t). They have been proposed in a different context in Busse & Clever ( 2000| , to 
guarantee mass conservation, J dzU(z) = 0. Finite Px{t), Py{t) appear only in the DNS 
of complex patterns in 

For the following discussions, a compact symbolic representation of the equations for 
the fields 9, /, $ is useful: 


C ^Vix, z, t) = £V{x, z, t) + 7V[F + U, V] 


(A3) 


with X = {x,y) and the symbolic vector V = [ 6 *,/,$]^. The symbol N stands for the 
nonlinear terms which consist of quadratic forms in 0 , /, $ and U. 

As an example, we show the explicit expressions for the linear terms of 9 and /. This 
allows us to immediately identify the corresponding components of the linear operators 
C,£: 


d 

9 = -RA2f^ 

at 


A7^9- Rsmj{US{z)dx)9, 


— i 9 (V^A 2 / = V^A 2 / - cosy A 26 ' + s'm-fdxd:,9 - — sm'yRF[UQ]f 


Pr 


Pr 


(A4o) 

(A46) 


with A 2 = {dxx + dyy). The term F[U§] = [Uq{z)V^ — d‘^^UQ{z)]dxA2 originates from 
the contribution of the basic mean flow Uq in (2.2) to the velocity u in (2.1&). Note that 


9, f are not coupled to $ in (A4). 


In general, equations (A3) are solved with the boundary conditions 9{z = ±1/2) = 0 
and / = i 92 / = <!) = [/ = 0 at z = ±l /2 which derive from the no-slip boundary 
conditions v{z = ±1/2) = 0. These conditions are automatically satisfied by the use of 


Galerkin expansions with respect to 2 . As in Busse & Clever (1992) we use for 9 the 
ansatz: 

M 

9{x,z,t) Sm{z)'dm{x,t)] (z) = sin(TO 7 r(z 1 / 2 )), (A 5) 


since Sm[z = ±1/2) = 0. For 4) and the secondary mean flow U{z,t) in (A2) also sine 


functions are used, while / is expanded in terms of the Chandrasekhar functions Cm{z) 
( Chandrasekhar|[l961| ) with (^^(±1/2) = dzCm{z = ±1/2) = 0. 


A.l. Linear Stability Analysis of the basic state 
The primary convectio n inst ability of the basic state corresponds to exponentially grow¬ 
ing solutions in time oA A3) in the linear regime (TV = 0). We use the ansatz V{x,z,t) = 


z, R) in (A3) to arrive at the following linear eigenvalue problem for a: 


aC{q, dx)V{q, z; R) = CV{q, z; R) = [A{q, d^) ± RB{q, 


(A 6 ) 


where the operators C,L{q,dz) etc. in Fourier space derive from the corresponding ones 
in position space (see (A3)) carrying a hat symbol via the transformation dx —t ig. In 
this paper we make use of Galerkin expansions ((see (A5)) to handle the z dependence. 
Thus, for instance ( |A 4[ ) is transformed into an algebraic linear eigenvalue problem of 
dimension 2M in the Fourier-Galerkin space. 

If cTmaxiR, Pr,j, q) is the eigenvalue with the largest real part in (A 6 ), then rolls be¬ 
come unstable when (7rnax{R, Pr, 7 , q) crosses zero. The standard procedure to determine 
the neutral surface R = Ro{Pr, 7 ; q) through the condition <7max{R = Rq, Pr, 7 , q) = 0. 
The minimum of RQ{q,"f) with respect to q gives the critical wavevector qc and the 
critical Rayleigh number R^, = T?o( 9 c) as function of Pr,j. If the frequency Wc = 
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ImlamaxiRc, Qc)] = 0 , the bifurcation of the basic state is stationary otherwise it is 
oscillatory. This method works in the ranges 7 ^ 12^ and 7 > I?'' in figure where 
the rolls are stable for R < {1 + einst)Rc{l)- Between 12° < 7 < 17°, where Ro{q,j) is 
not unique, we analyze <Jmax(,R, Br, 7 , q) = 0 using scans of 7 , q at fixed R to determine 
the upward-bent threshold of the wave roll instability, which limits the stability of the 
longitudinal rolls. 

It turns out, that eigenvalues a with Re[(Tmax\ ^ 0 are obtained by reducing the 
eigenvalue problem (A 6 ) to its 0, / part (A4|. Thus, the use of our Galerkin-expansions 
(A 51 reduces Eq. (A 6 ) to an algebraic linear eigenvalue problem with 2M x 2M matrices, 
which is analyzed using standard linear algebra codes (LAPACK). 


A.2. Secondary of instabilities of roll solutions 
For intermediate Pr, considered in this paper, the primary bifurcation to rolls with 
wavevector qc at R = Rc i s stat ionary. To construct the evolving finite-amplitude solution 
V = Vr for R > Rc from (^A we use the Fourier ansatz: 


i=N/2 


Vr{x,z) = 




(A 7) 


k=-N/2 


With respect to z, we introduce an additional Galerkin expansion (see (A 5)) of the {N+1) 
Fourier coefficients Vr{kqc, z). Furthermore, the Galerkin expansion of the mean-flow 
U using sine functions leads to 2M additional equations. Thus, we arrive at a system 
of 3M{N -I- 1) -I- 2M coupled nonlinear algebraic equations for all Galerkin expansion 
coefficients. This system is solved by Newton-Raphson methods. 


The iteration process is started from the weakly-nonlinear roll solution of (A3) char¬ 
acterized in Fourier space by the ansatz Vumiiqc, z) = A{qc,R)Vr 

_ _^ qc- For e > 0 

sion with respect to the small parameter e 


Z) — ^Iqc^ ti)Vcnax{.qcT Z^ R)j ^ '^raax is 
max at q = qc- For e > 0, a systematic expan- 
determines the amplitude A of Vwni via the 


solution of the amplitude equation: 


(Traax{.qc,R,Pr,^)A-cA\A\^ =0, with Re [cTmax ] OC e . 


(A 8 ) 


For stationary primary ILG bifurcations at intermediate Pr, both Umax and the cubic 
coefficient c are always real and positive and the bifurcation is thus forward (supercritical) 
at onset, i.e. \A\^ oc e/c increases continuously beyond the threshold e = 0 . 

To examine the linear stability of the roll solutions W (A7), we linearize (A3) with 
respect to an infinitesimal perturbation 6Vr{x, z,t) of Vr- Then we arrive at a set of 
coupled linear equations for the components of SVr, which are solved by the standard 
Floquet ansatz: 

k=N/2 

6Vr{x,zA)=e^*e^^'-^ ^ e^’^'^'^-^SVrikqc, z), (A9) 

k=-N/2 


Thus, we arrive at a linear eigenvalue problem for the eigenvalues A{s',qc,R). This 
is mapped, as before, to a linear algebraic problem using an Galerkin expansion of 
6Vr{kqc, z). The eigenvalue Aq with the largest real part defines the growth rate Ao(s') = 
Re[Ao(s')] of the perturbation 6Vr- 

Given that Ao(s') assumes its maximum XmaxiR, Pr,j) at s' = we then deter¬ 

mine the smallest Rayleigh number R = Rmstip, Pi") at which Xmax crosses zero for every 
R. In other words at P = Rinat the secondary instability of the rolls with wavevector 
qc occurs for given parameters Pr,j. When lOinst = Rinst)] 7 ^ 0 the secondary 
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7 M 

qc 

Rc 

7 

M 

qc 

Rc 

7 

M 

qc 

Rc 

0° 6 

3.1159 

1707.985 

90° 

6 

2.7920 

8504.200 

100° 

6 

2.7723 

9150.774 

8 

3.1162 

1707.824 


8 

2.8076 

8470.286 


8 

2.7985 

9108.171 

10 

3.1163 

1707.784 


10 

2.8059 

8476.495 


10 

2.7877 

9115.841 

12 

3.1163 

1707.771 


12 

2.8056 

8477.690 


12 

2.7871 

9117.068 


Table 1: Dependence of threshold values Qc and Rc for Pr = 1.07 of longitudinal rolls 
at 7 = 0° (horizontal cell), transverse rolls at 7 = 90° (vertical cell) and at 7 = 100° 
(heating from above) upon different Galerkin trnncation parameters M. 


bifurcation is called oscillatory, otherwise stationary. For an appropriate interpretation of 
the instability, we have to determine for a given the index kmax corresponding to the 
largest modulus \5Vr{kmax(lc^ z))\ of the expansion coefficients in (A9). This yields then 
the wavevector(s) Qinst of the dominant destabilizing mode(s) as Qinst = kmaxQ + ^'max- 
It turns out that kmax is always governed by the temperature component of 5Vr and 
that \kmax\ ^ 1 holds. 

The stability of roll patterns along the method described above has been intensively 
employed for the investigation of the standard isotropic RBC (7 = 0) by Busse and 
coworkers also for q ^ qg. They have constructed the Busse balloon (Busse & Clever 


1979 19961, which is the stability diagram of rolls with wavevector q in the R, Pr— 


parameter space. In our anisotropic ILC system the calculation of the full Busse balloon 
is possible if we replace qc by arbitrary q everywhere in equations (A 7 - A 91. 

Here we consider only the special case q = qc, where only finite s' perturbations 
have been found to be relevant. This has been discussed in detail in In particular in 
several cases the same maximal value of \5'V{kq, z))| is assumed at two different integers 
^ 2,^3 and thus two different dominant destabilizing modes q^st = Q 2 ,Q 3 exist. The 
resulting patterns are in addition often characterized by wavevector resonances of the 
form qi + 92 + 93 = 0 between the wavevector 91 = qc of the basic roll pattern and 
the 92 , 93 , which characterize the 3D patterns for R > Rinst{(lc,l) in Fourier space. The 
stability analysis yields also the relative phases of the three Fourier amplitudes. Because 
of translational invariance of the system in the x— and the y— directions two phases can 
be chosen to be zero without loss of generality; the third one is then determined by the 
ratio of the components of 514 (^ 29 o, z) and SVrik^qc, z) with the largest moduli. 

A.3. Accuracy of the stability limits 

The accuracy of all thresholds and results presented in this paper depend on the choice 
of the truncation parameter M of the Galerkin expansion with respect to z (see ( |A 5| )) 
and the truncation parameter N in Fourier space (see (A7, A9)). In this paper we have 
always chosen M = 8 and N = 5. By systematically increasing these parameters (see 
below) we have tested that this choice is sufficient to guarantee that the relative errors 
of all results given in this paper are below 0.1%. Even reducing the values to M = 6 and 
TV = 3 does not practically change the curves shown in this paper. 

In the table we give representative examples for linear threshold values Rc and Qc 
and their dependence upon increasing values for the Galerkin truncation parameter M. 
In addition, the determination of the codimension 2 point varies as 7 c 2 = 77.7857°, 
77.7462°, 77.7544° and 77.7560° when varying M as M = 6 , 8 ,10,12 respectively. 


Analogous convergence checks have been performed for the secondary instabilities of 
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the convection rolls. In general the data are more sensitive against changes of M than of 
N. We did all calculations with M = 8, which guarantees the same accuracy of the data 
as in the linear regime above. 

To guarantee an relative accuracy of better than 0.1% IV = 3 is sufficient at small e; for 
e = 0(1) one needs N ^ A. This conclusion is supported by the following representative 
data. The secondary instability towards wavy rolls in figure at emst = 0.8 occurs 
at 7 = 16.50° for iV = 5 and at 7 = 16.5288°, for iV = 4. For iV = 8 The LSO 
instability at 7 = 17° is characterized by einst = 1.04423 with the circular oscillation 
frequency ujinst = 10.20808 and Floquet vector as Sx = 1.27929 and Sy = gcol"^- With 
the smaller = 4 we find only small changes with emst = 1.0436, = 10.2038 and 

Sx = 1.278990, Sy = qco/‘^- Finally we mention the knot instability of transverse rolls at 
7 = 83°. For N = 5 we find emst = 0.07134 with Sy = 3.1179, which remain unchanged 
for = 4. 


Appendix B. Direct simulations of the OBE in ILC 

Direct simulations of the OBE in (2.1) are in general confined to a rectangle in the x — y 
plane with the lateral extensions Lx,Ly using periodic boundary condition V{x,y,z) = 
V{x + Lx, y + Ly, z). Thus, we transform to Fourier space by introducing a discrete 2D 
Fourier transformation oi V on a Nf x Nf grid with mesh sizes Ag^,, Aqy in the x- and 
y-directions: 

V{x,z,t) = {q, z,t) where q = {{kAqx,lAqy)} with —Nf/2 < {k,l) ^ Nf/2, 

(Bl) 

and /Aqx = 2Tr/Lx,Aqy — 2 tt jLy. Reality of V{x,z,t) implies the condition V{q) = 
V {—q)*- With respect to z, we use Galerkin expa nsions with the truncation parameter M 
as before. The quadratic nonlinearities A in (A 31 are treated by standard pseudospectral 
methods (see e.g. Boyd (2001)). Substituting the Fourier ansatz (Bl) into (A3) and 


projecting on the respective Galerkin modes one arrives at a system of 3M x Nj coupled 
ordinary differential equations for the evolution of all the combined Fourier-Galerkin 
expansion coefficients. In addition, (A2) is mapped into a system of 2 M equations 
for the Galerkin coefficients of the secondary meanflow U. Semi-implicit time stepping 
methods, as sketched in the following subsection, are used to compute the time evolution 
of all our fields. 

For all DNS shown in section we have used Lx = Ly = n^Xc with the critical 
wavelength Ac = 2,'Klqc and the truncation parameters Ay = 256, M = 8 and ul = 20. 
Note that the number of roll pairs (black and white stripes) of the underlying 2D roll 
patterns directly reflects n^. 

As also evident from the previous section our Fourier coefficients decay quickly with 
increasing |q|. The truncation parameter Ay = 256 in the simulations corresponding to 
A = (124/20) > 5, is used for the stability analysis of rolls in Appendix A.2 Thus it 


is not surprising that the results of the Floquet analysis are reproduced in the DNS. 
Increasing Ay and/or decreasing ul corresponds to keeping Fourier modes with larger 
|q|. We have checked that all the typical scenarios discussed in )|^are recovered. 

The goal of most simulations i n (j^ was to validate the secondary instabilities of rolls 
originally obtained on the basis of |A(2[ which lead to strictly-periodic 3D patterns. Thus, 
we have performed the DNS on minimal domains in the plane, where one side length 
was given as Ac while the other was determined by the wavevectors of the dominant 
destabilizing modes q 2 , <73 introduced in Appendix|A.2[ The data have then been mapped 
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to a larger domain in the x — y plane by periodically extending the minimal domains for 
visualization. 


B.l. Exponential Time Differencing method 


Our starting point is (A3), where the components of V(q,z,t) in (Bl) are expanded 
into the appropriate Galerkin modes like in (A5). In the resulting Fourier-Galerkin (A3) 
can be written as: 


V{t)=AV(t)-N with N =C-^N, A = C-^C, 
at 


(B2) 


since the matrix C is not singular. Equation (B 2) allows for the formal solution: 

N{t')dt'. (B3) 


/ I 


Approximating N{t') by the leading terms of the Taylor expansion about the lower 

t' -\-t and 


limit t of the integral in (B 3) followed by the variable transformations t' 

N{t) - N{t-dt) 


subsequently r' = rdt one arrives at: 


V{t + dt)=e^‘^^V{t)-dt e 


{l — r)Adt 


N(t) T dt ■ 


dt 


dr. (B 4) 


In Koikari (2009), one finds for an arbitrary matrix M, the following dehnition of the 

f g(l-r)M k=l^2.... (B5) 

• Jo 


matrix functions (^^(M): 


1 

(^ 1 )! 


Thus, (B4) can be rewritten as: 


V{t-\-dt)=(j)Q{Adt)V{t)—dt(j)i{Adt)N{t) — dt^(l)2{Adt) | (B6) 


Note, that both V and the secondary meanflow U (see A2)) are calculated using the 
time exponential method. 

The time step ping scheme described in (B 6) incorporates the convergence to stationary 
solutions Vs of (B2| which have to fulfill AVg — Ng = 0. This can be proven using the 
recurrence identities of the matrix operators ^^(M) as 

MM) = - I] ■ (B7) 

It should be remarked that the matrix exponentials could be also treated by using 
a spectral representation of A (B2) in terms of its direct and adjoint eigenfunctions. 
In this way, one makes immediate contact to the method used in Pesch (1996). While 
this procedure has been successfully applied in a series of papers on complex patterns 
in standard RBC (Bodenschatz et al. 2000 Egolf et al. 2000), its application to ILC 
requires particular care and is less robust, since the spectral properties of the operator 
A are complicated ( Rudako^|1967 Ghen &: Pearlstein]|1989 ). 


Appendix C. Additional remarks on the linear stability calculations 

In (i the linear instability of the basic state, against rolls with wavevector q = 
q{cosip,sui'ip) was discussed. This case corresponds to the existence of an eigenvalue 
a of (A6| with Re\a\ > 0]. In the following, we will exclusively concentrate on the 6,f 
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components obtained from Fourier transformation of (A4)), since < 0 for eigenval¬ 
ues of the separated $— equation. Since (A4) are invariant against the transformations 
[x,z) —> —{x^z) and separately (y —>■ —y), it is sufficient to restrict "0 to the interval 
(0 ^ ^ < 90°). As already mentioned in and documented in figure]^ for Pr = 1.07, 
it is sufficient to investigate only the special cases either ip = 90° (longitudinal rolls) or 
Ip = 0° (transverse rolls). A proof can be found in Gershuni & Zhukhovitzkii ( 1969| ). In 
subsection |C.2[ we will present our own very short version. 

C.l. Linear stability results for small and large Pr 
We have reproduced some of the earlier results in the literature as validation of our nu¬ 
merical methods. In general, the codimension 2 point 7 = 7 c 2 , where the critical Rayleigh 
numbers i?c( 7 ) of longitudinal rolls and i?c( 7 ) of the transverse ones are equal, moves 
continuously towards 7 = 0 for decreasing Pr . Below Pr < 0.264 (more precisely cal¬ 
culated in Fujimura & Kelly (1992)), only the primary bifurcation to transverse rolls 


prevails even at incremental inclinations from the flat 7 = 0 state. 

On the other hand, for Pr > 1.07, the codimension 2 point moves continuously to¬ 
wards 7 = 90® and eventually the bifurcation to transverse rolls is relevant only in the 
range 90° < 7 < 180°. In addition, the primary bifurcation to transverse rolls becomes 
oscillatory at large Pr. For the vertical case (7 = 90°) we calculate that this happens 


at Pr > 12.45, in agreement with Fujimura & Kelly (1992). The oscillatory bifurcation 


has been perfectly reproduced by our own calculations. At a slightly larger Pr = 12.7 
again at 7 = 90®, we find Rc = 88220, Wc = 504.6 which is in excellent agreement with 


Bergholz (1977). 


C.2. Bifurcation of oblique rolls 


It is useful to exploit a certain form invariance of the linear equations (A4a A46) in 
Fourier space using the transformation / —>• //i?. It turns out that the Rayleigh number 
R, the inclination angle 7 and the oblique roll angle ip appear first in the combination 


i?sin 7 cos■)/' due to the contributions of C/q (see ( 2 . 2 )) and furthermore in (A46) as 


Rcosj. In addition, these equations depend only on |q 

Thus, the critical Rayleigh number R}f.{l) of the longitudinal rolls {ip = 90°) is deter¬ 
mined by Rcil) cos 7 = i?c (7 = 0) = 1707.824. At the codimension 2 point 7 c 2 , we have 
R’‘c{ 1 c 2 ) = RI{ 1 c 2 ), where R^{j) denotes the critical Rayleigh number of the transverse 
rolls. Thus, 7 c 2 is determined by the transcendental equation i?‘( 7 c 2 ) cos( 7 c 2 ) = Rco 
which has always a unique solution 0 < 7 c 2 < 90° for 0.246 < Pr < 12.45. 

The form invariance of the 6, f equations implies a direct relation between transverse 
eigensolutions {ip = 0) of (A 6 ) for given values of 7 90°, \q\,a,Pr and the oblique 

eigensolutions {ip ^ 90°) for the same values of \q\,a,Pr but at a different inclination 
angle F ^ 7 and for a Rayleigh number i? = i?‘(F) 7 ^ R°^. In detail, we have: 


R°° sin 7 cos Ip = R*{T)smT, 


R°^ cos 7 = cos F 


tan( 7 ) cos{ip) = tan(F). (C 1) 


In the special case F = 7 = 90° equation (Cl) simplifies to R°° cos ip = R^ .In Gershuni 


& Zhukhovitzkii (1969) the authors have arrived to analogous relations. As a general 


consequence an explicit analysis of linear oblique rolls is not necessary, since they can be 


determined from the transverse rolls on the basis of (Cl). 


Equations (Cl) are useful to characterize the stationary bifurcation to oblique rolls 
(cr = 0) at medium Pr, which is characterized by the curve R = R°^{'y,ijj) as function 
of the inclination angle 7 . This curve crosses the longitudinal threshold curve Rc{'j) 
at the codimension 2 point 7 °|('!/') < 90°, which is thus determined by the relation 
^f(7°2,V’) = i?co/cos(7,1)- 
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Thus, (Cl) predicts that there exists a certain angle T' which allows for expressing 
by RliT') as follows: 

i’) cos(7 °2) = (r') cos(r') = tan(7;;'2) cos(V’) = tan(r'). (C 2 ) 


From the first equation in (C 2), we conclude F' = 7^2 and from the second one > 7 c 2 
for {0 < il) < 90°), such that for 7 > 7 c 2 pure transverse rolls prevail at onset. From 
( |C 2[ ), is also obvious that the dip in the gc( 7 )—curve in the transverse case {-ijj = 0) is 
mapped to corresponding ones for ip ^ 0 with the same \qc\] this immediately explains 
their equal heights in figure 

Though not relevant for the onset of convection, the sudden, strong increase of i?*{ 7 ) 
curve at 7 « 26.5° in figure]^ is interesting. For Pr > 1.75 even discontinuous jumps 
in i?* and q* at 7 « 30° develop. These are associated with the phenomenon of two 


disconnected neutral curves as seen in figure 6 in Chen & Pearlstein (1989)). 
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